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ABSTRACT 

We report the discovery of 31.18 ms pulsations from the INTEGRAL source IGR J14003-6326 using the 
Rossi X-ray Timing Explorer (RXTE). This pulsar is most likely associated with the bright Chandra X-ray point 
source lying at the center of G3 10.6-1.6, a previously unrecognised Galactic composite supernova rem- 
nant with a bright central non-thermal radio and X-ray nebula, taken to be the pulsar wind nebula (PWN). 
PSR J 1400-6325 is amongst the most energetic rotation-powered pulsars in the Galaxy, with a spin-down lu- 
minosity of E = 5.1 x 10 37 erg s _1 . In the rotating dipole model, the surface dipole magnetic field strength is 
B s = 1.1 x 10 12 G and the characteristic age t c = P/2P = 12.7 kyr. The high spin-down power is consistent 
with the hard spectral indices of the pulsar and the nebula of 1.22 ± 0.15 and 1.83 + 0.08, respectively, and a 
2-10keV flux ratio Fpwn/Fpsr ~ 8. Follow-up Parkes observations resulted in the detection of radio emission 
at 10 and 20 cm from PSR J1400-6325 at a dispersion measure of ~ 560 cm -3 pc, which implies a relatively 
large distance of 10 ± 3 kpc. However, the resulting location off the Galactic Plane of ~ 280 pc would be 
much larger than the typical thickness of the molecular disk, and we argue that G3 10. 6- 1.6 lies at a distance 
of ~ 7 kpc. There is no gamma-ray counterpart to the nebula or pulsar in the Fermi data published so far. A 
multi-wavelength study of this new composite supernova remnant, from radio to very-high energy gamma-rays, 
suggests a young (< 10 3 yr) system, formed by a sub-energetic (< 10 50 ergs), low ejecta mass (M e j ~ 3 M ) 
SN explosion that occurred in a low-density environment (no ~ 0.01 cm -3 ). 

Subject headings: pulsars: individual (PSR J1400-6325) — X-rays: individual (IGR J14003-6326, 
G3 10. 6- 1.6) — supernova remnants — gamma rays: observations 



1. INTRODUCTION 

The number of known supernova remnants (SNRs) in the 
Galaxy has increased significantly over the last few years, 
mainly due to a new generation of radio and X-ray instru- 
ments of unprecedented sensitivities and angular resolutions. 
In particular, these observations, com bined with Galactic 
Plane surve ys (e. e. iBrogan et al. 112006 ), and targeted obser- 
vations (e.g. iGelfand & Gaenslerl2007l:lGaensler et al. 120081) 
have increased the fraction of SNRs found to harbor an ener- 
getic pulsar powering a wind nebula (the so-called composite 
SNRs). Moreover, deep observations toward compact PWNe 
have proven to be succes s ful in detecting their powe ring pul- 
sars (ICamilo et al. 1120091: iGotthelf & Halpernll2009h . provid- 
ing constraints on their energetics, spin evolution, and birth 
parameters. Along with the current generation of high/very- 
high energy (VHE; > 100 GeV) instruments, these observa- 
tions are of prime importance for understanding the struc- 
ture and evolution of these sources, and the underlying ac- 
celeration mechanisms which occur close to the pulsar, and at 
the relativistic and non-r elativistic shock fronts b ounding the 
PWN and the host SNR dGaensler & Slane 20(j3). 
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The soft y-ray source IGR J 14003-6326 was discov- 
ered in a deep INTEGRALfTBlSflSGRl mosaic of the Circi- 
nus region as a persistent so urce at the mCrab level 
dKeek. Kuiper & Hermsenl 120061) . A Swift /XRT survey of 
INTEGRAL sources located IGR J14003-6326 to 4" in 2- 
10 k eV X-rays, but no c onclusion was reached on its na- 
ture dMalizia et al. 1 12007). A follow-up Chandra survey of 
unidentified IGR sources reported a P WN within a ~ 3' di - 
ameter nearly circular emission nebula dTomsick et al. 120091) . 
These authors presented the 0.3-10 keV spectrum of the to- 
tal emission by an absorbed power-law with a relatively hard 
photon index (F = 1.82 + 0.13) and a large column density 
N H ~ 3 x 10 22 cm -2 . 

In this paper, we report the discovery with the Rossi 
X-ray Timing Explorer (RXTE) of 31.18 ms pulsations to- 
ward IGR J14003-6326, and the radio detection in follow- 
up observations using the Parkes telescope. We present 
a multi-wavelength study of this new Galactic composite 
SNR, G310.6-1.6, using radio, X-ray, and gamma-ray data. 
PSR J1400-6325 is one of the most energetic in the Galaxy 
and powers a wind nebula whose broadband non-thermal 
synchrotron emission is measured in radio and X-rays. We 
present a spatial and spectral analysis that shows evidence for 
a young SNR, formed from a sub-energetic, low ejecta mass 
SN explosion that occurred in a low-density environment. We 
also discuss the implications of the lack of high-energy (HE) 
gamma-ray emission in the Fermi data. 

2. OBSERVATIONS AND RESULTS 
2.1. Chandra 

We have performed a reanalysis of the Chan dra data of the 
IGR J14003-6326 field dTomsick et al. Il2009l) taken on 2008 
June 29 UT, using the Advanced CCD Imaging Spectrometer 
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Fig. 1. — Left: Chandra/ ACIS-I count rate image of G310.6-1.6 in the 0.5-10keV energy r ange, smoothed with a gaussian ofcr = 5" (linear scale). The dashed 
ellipse shows the intrinsic size of the radio MGPS-2 extended source cataloged at 843 MHz I Murphy et al. 2007, see section l2~5l . The bright, central, point-like 
source, is better revealed in the inset unsmoothed image of the central regions of the SNR. Right: Surface brightness profile of the Chandra/ACIS emission, 
centered on the point-like source. The background level is depicted by the horizontal solid line. The dashed lines delineate the regions of interest, namely the 
PSR, the PWN, the putative X-ray halo, and the host SNR. 



(ACIS. iGarmire et al. ||2003|) and operating in TIME/FAINT 
exposure mode. The livetime amounts to 5077 s, corrected 
from the 1.3% readout dead time. The maximum pixel count 
rate within the source is 0.015 ct s _1 , so pile-up effect could 
be safely ignored. All the data reduction and analysis were 
performed with the Chandra Interactive Analysis of Obser- 
vation software (CIAO) version 4.1.2, using the CALibra- 
tion DataBase (CALDB) v4.1.3. Figure Q] (left) shows the 
0.5-10 keV Chandra image centered on IGR J14003-6326, 
which exhibits a characteristic morphology of a composite 
SNR, as revealed by the radial profile of the surface bright- 
ness (Figure Q] right). In particular, the excess emission at 
1.1-1.2' from the center is found to be the clear signature of 
the host shell-type SNR. We have considered three regions of 
interest, namely the bright point-like source (PSR), the dif- 
fuse surrounding emission (PWN), and the shell-like struc- 
ture (SNR). As shown on the profile, some faint emission 
seems to extend up to the shell-type structure, at distances 
of ~ 1' from the PSR. In view of the large Nh, this could be 
partly explained by a dust-scattering halo, s imilar to those ob- 
served in other PWNe such as G21. 5-0.9 dSlane et al. H2Q00t 
lSafi-Harbetai~ll200Tt iBocchino et al. 1120051) . Unfortunately, 
the lack of statistics prevents us from drawing any firm con- 
clusion about the nature of this extended X-ray emission. The 
source-free region used to extract the background was taken 
on the same CCD chip as the source, in a circle of 1.45' in ra- 
dius. Apart from the PWN emission, the low count levels led 
us to fit the three components simultaneously, with a common 
hydrogen column density Nh = ( 2.09 + 0. 12) x 10 22 crrT 2 . As 
noticed bv lTomsick et al. I (120091) . this value is of the same or- 
der as the Galactic valufl and suggests a large distance and/or 
local absorption. 

The best-fit position of the point-like source 
is R.A.(J200 0) = 14 h 00 m 45.69 s , Decl.(J2000) = 
-6325'42.6" (ITomsick et al. I 12009b . For the spectral 

7 We here make use of t he Anders & Grevesse 1 1989) abundances, rathe r 
those of IWilms. Allen & M cCrav I 20Q3), as done in ITomsick eTaTl <2"0Tj9T) . 
Using the latter abundances yields compatible Nh ~ 3 X 10 22 cirT 2 . 



analysis, we chose an integration region of 1.5" in radius 
centered on the source and estimated the source background 
using counts extracted from a concentric annulus of radii 
1.5" and 2". Fits to the background subtracted spectrum 
are not tightly constrained, even though a power-law gives 
slightly better results than any thermal model. The spectrum 
is hard, with a best-fit power-law index of 1.22 + 0.15 
and an unabsorbed 2-10 keV flux of (1.95 ± 0.5) x 10~ 12 
ergs cirT 2 s _1 . The PWN integration region was taken in 
an annulus centered on the PSR position, with an inner 
radius of 1.5" and an outer radius of 0.65' (see Fig. Q] 
right). The best-fit power-law index and the unabsorbed 
2-10 keV flux are found to be 1.83 + 0.09 and (1.51 ± 
0.20) x 10~ 11 ergs cirT 2 s -1 , respectively. Therefore, the 
PWN is by far the dominant component in X-rays and these 
sp ectral results are in clo se agreement with those obtained 
by ITomsick etaD (120091) on the total emission. The SNR 
integration region was defined as an annulus, with inner and 
outer radii of 0.98 and 1.27', respectively. As for the PSR, 
due to the low statistics, the SNR spectrum is rather poorly 
constrained, and does not allow us to clearly distinguish 
between thermal and power-law models, although no strong 
X-ray lines are present. A power-law fit gives a rather soft 
spectrum, with a best-fit spectral index of 2.56 ±0.18 and an 
unabsorbed 2-10keV flux of (1.0 + 0.2) x 10~ 12 ergs cm" 2 
s . Whatever is the nature of the X-ray emission, the flux 
measured with Chandra can be considered as an upper limit 
on the synchrotron emission. We then converted the best-fit 
power-law spectrum to a surface brightness at 4 ke V of ~ 5 x 
10 32 ergs pc~ 2 s sr _1 , as this will be compared to theoretical 
predictions in section[3] 

2.2. RXTE 

The X-ray flux measured from the Chandra point source 
in G3 10. 6- 1.6 is sufficient to search for the expected pulsa- 
tions using RXTE. A pair of 44 ks exposures were obtained 
a year apart on 2008 September 29 UT and 2009 September 
30. Data were collected with the Proportional Counter Ar- 
ray (PCA; iJahoda et al. II1996T) in the GoodXenon mode with 



Discovery of PSR J1400-6325 associated with IGR J14003-6326 



3 



an average of 1.5 and 2.2 out of the five proportional counter 
units (PCUs) active, for the two observations, repectively. In 
this mode, photons are time-tagged to 0.9 fis and have an ab- 
solute uncertainty better than 100 fis. The effective area of 
five combined units is 6500 cm 2 at 10 keV with a roughly 
circular field-of-view of ~ 1 FWHM. Spectral information is 
available in the 2-60 keV energy band with a resolution of ~ 
16% at 6keV. The standard time filters were applied to the 
PCA data, which rejects intervals of South Atlantic Anomaly 
passages, Earth occultations, and other periods of high par- 
ticle activity. The photon arrival times were transformed to 
the Solar-system barycenter in Barycentric Dynamical Time 
(TDB) using the JPL DE200 ephemeris and the Chandra co- 
ordinates for the point source. 

2.2.1. Timing analysis 

We restricted the RXTE timing analysis to photon data with 
energies in the 2-20 keV range (PCA channels 2-50) recorded 
in the top Xenon layer of each PCU, to optimize the signal- 
to-noise ratio. For the first observation, we further excluded 
one 4 ks non-contiguous data segment at the beginning of the 
run, resulting in 33.4 ks of filtered time spanning 45 ks. A 
2 26 Fast Fourier Transform (FFT) of the data binned in 1 ms 
steps revealed a significant signal of power S = 55 at a period 
of P = 31.1 ms. For the 6.7 x 10 7 search elements this corre- 
sponds to a false detection probability for a blind search down 
to 1 ms of p = 7.6 x 10~ 5 (99.9924 % Confidence Level). We 
refined this signal using a Z 2 statistics which yields a signal 

power Z 2 = 53.76 at a period of P = 3 1.1 80373(2) ms. The 
lcr uncertainty on the last digit is given in parentheses. Fig- 
ure |2] displays the pulse profile in the 2-20 keV band folded 
at this period; the profile is narrower than a sinusoid, some- 
what boxy, with a well deliniated off-pulse region. We do 
not see any evidence for energy dependence of the pulse pro- 
file when subdividing the 2-20 keV band. The signal does 
not match any cataloged pulsar in the field-of-view, in par- 
ticular the pulsars PSR J1403-6310 (P = 0.39910 s) and 
PSR J1401-6357 (P = 0.84279 s). In order to ensure that 
we have not picked up an aliased signal from these PSRs, the 
data was rebarycentered at their respective coordinates and 
searched for the known pulse period. No signal was detected 
apart from the original one reported above. 
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PULSE PHASE 

Fig. 2.— PSR J1400-6325 folded light curve in the RXTE PCA 2-20 keV 
band. Two cycles are shown for clarity. Phase zero is arbitrary. 



TABLE 1 

Measured and derived parameters for PSR J 1400-6325 



Parameter Value 



R.A. (J2000) 14 h 00 m 45.69 s 

Decl. (J2000) -63°25'42.6" 

Period, P, (ms) @ MJD 54738 31.180373(2) 

Period, P 2 (ms) @ MJD 55105 31.1816011(3) 

Period derivative, P 3.890(6) X 10" 

Dispersion measure (cm" 3 pc) 563 ± 4 

Flux density at 10 cm (/iJy) ~ 110 

Flux density at 20 cm (/ijy) ~ 250 

Flux density at 50 cm (//Jy) < 500 

Characteristic age, t c (kyr) 12.7 

Spin-down luminosity, E (ergss" 1 ) . 5.1 X 10 37 

Surface dipole magnetic field, B s (G) 1.1 X 10 12 

DM-based distance (kpc) 10 ± 3 



Note. — lcr uncertainties given. 

Based on a probable detection of an highly energetic PSR, 
we requested a second RXTE observation to confirm this re- 
sult and determine the PSR spin-down rate. For the 2009 data, 
we followed the same proceedures outlined above and con- 
structed a Z 2 periodogram around the expected signal. This 
yielded a significant signal at P = 31 .181601 1(3) ms with 
power Z 2 = 116.28, corresponding to a negligible proba- 
bility of false detection. This result confirms the previous 
detection and provides a period derivative measurement of 
P = (3.890 ± 0.006) x 10~ 14 . The one sigma error is de- 
rived from the root sum of the variance of the period uncer- 
tainties. The ephemeris is presented in Table[TJalong with the 
inferred spin-down parameters. The derived spin-down power 
£ = 4tt 2 IP/P 3 = 5.1 x 10 37 h 5 ergs s" 1 (7 45 is the moment of 
inertia in units of 10 45 g cm 2 ) makes PSR J1400-6325 one 
of most energetic Galactic pulsars known, a close tie with 
several recently reported energetic rotation-po wered pulsars 
dCamilo et al"1l2009t iGotthelf & Halpernll2009l) . In the dipole 
pulsar model, the characteristic age t c = P/2P = 12.7 kyr, 
assuming P » Pq (the initial PSR period), and the surface 
dipole magnetic field B s = 1.1 x 10 12 G. Based on the as- 
sociated PWN and SNR, its estimated luminosity, and its 
spin-down properties, PSR J 1400-6325 is clearly a Crab-like, 
young energetic rotation-powered PSR responsible for the ob- 
served PWN. 

The pulsed emission of PSR J 1400-6325 is well isolated 
in phase with nearly constant emission in-between. This al- 
lowed us to take advantage of phase-resolved spectroscopy to 
obtain a spectrum of the pulsed flux using the off-pulse emis- 
sion as a near perfect representation of the background plus 
unpulsed spectrum. We generated phase-resolved spectra in 
20 phase bins using the RXTE WOOL fasebin, aligning the 
two observations to a common phase and summing the result. 
From this combined spectral file, we generated on- and off- 
pulsed spectra as input into the XSPEC spectral fitting pack- 
age. The response matrix for the phase-resolved spectra was 
generated using the same photon file, selecting only photons 
from the top layer of each PCU (LR1) to match the fasebin 
output. We divided the data in to two phases, 0.5 cycles apart 
to characterise the on-pulse signal and the off-pulse signal, re- 
spectively. We fitted the on-pulse spectrum with an absorbed 
power-law spectral model over the 2-10 keV band, outside of 
which the background dominated the pulsed signal. The col- 
umn density was fixed to the Chandra value as RXTE, with its 
higher energy band, is unable to constain this parameter. For 
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an N H = 2.09 xlO 22 cirT 2 the best-fit photon index is 2.0^1 
with ax 2 = 1.2 for 17 degrees-of-freedom (see Fig. [3}. The 
measured 2-10 keV pulsed flux is 3.0 x 10~ 13 ergs cirT 2 s _1 
(absorbed), and represents ~ 18% of the total flux from the 
Chandra source. 




(KEV) 



Fig. 3.— PSR J1400-6325 pulsed spectrum in the RXTE PCA 2-20 keV 
band fitted to an absorbed power-law model given in the text. The residuals 
from the best fit model are show. The off-peak spectrum is used as back- 
ground (see text). 



2.2.2. Spectral analysis 

In order to extract the total RXTE/PCA spectrum, domi- 
nated by the PWN, we also accumulated spectra restricted to 
the top layer of PCU 2. The background was estimated using 
the faint model. Given the sky position of IGR J14003-6326, 
X-ray emission from the Galactic ridge can contri bute to the 
PCA spectrum. Following p revious works (e.g. iPrat et al71 
12001 iRodriguez et"al~lr2 009), we estimated its level accord- 
ing to the lValinia & Marshalll (11998b model, and by assuming 
that Chandra and RXTE/PCA are perfectly cross-calibrated. 
After removing the Galactic X-ray contribution, we found that 
a power-law fit to the 3-20keV PCA spectrum gives a Nh sim- 
ilar to that found for the three Chandra components, though 
not well constrained, (2.25 + 1.1) x 10 22 cm~ 2 . Freezing it 
to the value measured with Chandra leads to a spectral in- 
dex of 2.19 ± 0.10 and an unabsorbed 3-20 keV flux of (1.57 
+ 0.3) x 10~" ergs cirT 2 s -1 . The spectrum measured with 
PCA turns out to be softer than that measured with Chan- 
dra/ AOS in the 2-10 keV band. This can be explained by 
the presence of a break or steepening at high energies in the 
PWN spectrum. Thus, we also fitted the PCA spectrum with 
a broken power-law, with Nh and a photon index before the 
break frozen to the Chandra values, and found a marginal ev- 
idence for such a break at ~ 3 keV, with a photon index at 
higher energies of 2.21 ± 0.09, though with a slightly larger 
X 1 . This emission detected by RXTE/PCA is most likely re- 
lated to G310.6-1.6, and mainly to the PWN diffuse emis- 
sion, since no other significant X-ray source lies within the 
PCA field-of-view, at 2 from IGR J 14003-6326. 

2.3. INTEGRAL IBIS/ISGRI 

We analyzed all the public IBIS/ISGRI data, from revo- 
lution 46 to 713, at less than 12 from IGR J14003-6326, 
with the Off-Line Scientific Analysis (OSA) software 



dGoldwurm et al. 112003b . version 7.0. After removing all the 
noisy pointings, as defined in iRenaud et al. I (120061) . the to- 
tal effective time amounts to ~ 3.5 Ms. Mosaic images were 
produced in six energy bands from 18 to 65 keV with an ad- 
ditional band between 82 and 150keV, in order to extract the 
non-thermal continuum source spectrum. Furthermore, the 
very small apparent size of G310.6-1.6 (O ~ 2.3' = 4.7 dj pc, 
with d = 7 dj kpc) is reminiscent of a young SNR. To search 
for evidence of 44 Ti line emission at 67.9 and 78.4 keV, ex- 
pected from a young remnant, we also generated mosaic im- 
ages in the 65-71, 71-75 and 75-82 keV energy bands. These 
lines emerge from the radioactive decay of the short-lived ^Ti 
nucleus (with a weighted-average lifetime of 85 + 0.4 yr, see 
Ahmad et al. 2006, and references therein). Alt hough firmly 
detec t ed in only the - 330 yr old Cas A SN R (llvudin et al. I 
U994tlVink et al. 11200 UlRenaud et al. II2006I) . this radioactive 
element can provide invaluable constraints on the key param- 
eters of the SN explosion. 

In the 18-50 keV mosaic image, the best-fit position of 
IGR J 14003 -6326 is found to be consistent with the results 
obtained by iBird et aiTl (12007b . and compatible with a point- 
like source concident with G3 10. 6- 1.6. The continuum emis- 
sion is safely detected with IBIS/ISGRI up to 65 keV and the 
spectrum is fit by a power-law with a photon index of 2.33 
± 0.29 and a 20-100 keV flux of (1.24 + 0.17) x 10 -11 ergs 
cirT 2 s _1 . No evidence of 44 Ti line emission was detected in 
the direction of IGR J 14003-6326, and we set a 3 <x upper 
limit of 1.5 x 10~ 5 ph cm~ 2 s _1 , after combining the two ^Ti 
energy bands. 

2.4. Spectral properties of G310.6—1.6 
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Fig. 4. — Spectra of the different components of G3 10. 6- 1.6 measured with 
Chandra (in black). The SNR spectrum was divided by 2 for sake of clar- 
ity. The RXTE/PCA and ZW£Gft4L/IBIS/ISGRI spectra are shown in blue 
and red, respectively. Solid lines represent the best-fit models (see Table 0. 
An absorbed power-law was used to fit simultaneously the individual Chan- 
dra components, while an absorbed broken power-law was fit to the Chan- 
dra spectrum of the whole X-ray emission, together with the RXTE/PCA and 
IBIS/ISGRI spectra. Residuals are shown in the lower panel. 



Spectra of the different components measured with Chan- 
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TABLE 2 

Imaging-spectroscopy of G310.6-1.6, PSR J1400-6325 and its pulsar wind nebula 





Nn.' 


power- 

r 


■law 
Flux b 


x 2 /y 


Nk* 




broken 


power-law 

r 2 


Flux b 


2; 


PSR (pulsed) RXTE 


2.09 (-) 


2.0 (0.40) 


0.36 (0.2) 


20.4/17 














PSR (total) Chandra 


2.09(0.12) 


1.22 (0.15) 


1.95 (0.5) 


182.5/191 














SNR Chandra 


2.09(0.12) 


2.56 (0.18) 


1.00 (0.2) 


182.5/191 














PWN Chandra 


2.09(0.12) 


1.83 (0.09) 


15.1 (2.0) 


182.5/191 














RXTE/PCA 


2.09 (-) 


2.19 (0.10) 


15.7 (3.0) 


31.8/30 


2.09 (-) 


1.83 (-) 


3.0 (2.0) 


2.21 (0.09) 


15.7 (3.5) 


31.8/29 


IBIS/ISGRI 




2.33 (0.29) 


12.4(1.7) 


7.2/7 














Total d 


2.46 (0.08) 


2.16(0.03) 


13.4 (0.8) 


272.6/246 


2.13 (0.13) 


1.90 (0.11) 


5.6 (0.7) 


2.30 (0.07) 


10.6(1.5) 


257.1/244 


PWN C 










2.13 (0.14) 


1.90 (0.10) 


6.0 (0.5) 


2.59 (0.11) 


5.3 (0.8) 


257.1/244 



Note. — The uncertainties quoted are at the 68 % confidence level. 

In units of 10 22 citT 2 , assuming the Anders & Grevesse (1989) abundances 

2-10 {Chandra), 3-20 (RXTE/PCA), and 20-100 keV (IBIS/ISGRI, Total and PWN) unabsorbed fluxes, in units of 10~ 12 ergcnT 2 s" 



Break energy in units of keV 
Total spectrum in the 0. 8-100 keV energy range 

PWN spectrum in the 0. 8-100 keV energy range, after removing the PSR and SNR contributions 



dm, together with those of RXTE/PCA and IBIS/ISGRI are 
shown in Figure |4] and spectral model fits are reported in Ta- 
ble El In the 2-10 keV band, the PWN accounts to ~ 85 % 
of the whole emission. A simultaneous fit on the three spec- 
tra (Chandra total emission, RXTE/PCA and IBIS/ISGRI) in- 
dicates that a broken power-law is better at describing the 
0.8-100 keV total emission. According to the F-test, this 
model is favored over a pure power-law at the 3.5 cr confi- 
dence level (probability of 2 x 10~ 4 ). The break energy is 
measured at 5.6 + 0.7 keV, and the photon index at higher en- 
ergies nicely matches that measured by IBIS/ISGRI. We also 
performed a spectro-imaging analysis of the PWN, in order 
to search for any spectral softening at larger distances from 
the PSR, as expected from efficient synchrotron burn-off of 
high -energy electrons accelerated in the central regions (see 
e.g. ISlane et al. Il2000l in the case of G21. 5-0.9). We fixed 
the Nh to the best-fit Chandra value, and found a 2.5 cr evi- 
dence of a spectral softening, from r = 1.68 ± 0.09 at 5" from 
the PSR, to 2.09 + 0.09 at the edge of the nebula. In order to 
measure the 0. 8-100 keV PWN spectrum, we first constrained 
the PSR spectral index and 2-10 keV flux to be ~ 1.4 and 
1.7 x 10~ 12 ergs cirT 2 s _1 at most, in order not to violate the 
IBIS/ISGRI 3 cr upper limit in the 82-150keV energy range, 
while being still consistent with the Chandra values, at the 1 cr 
confidence level. Then, we performed a second simultaneous 
fit on the three spectra (Chandra total emission, RXTE/PCA 
and IBIS/ISGRI), by freezing the 0.8-10 keV PSR and SNR 
spectra to the best-fit values, and extracted the minimum 0.8- 
100 keV spectrum of the PWN, whose parameters are pre- 
sented in Table [2] (last row). We find that a broken power-law 
fit gives a break energy of 6.0 ± 0.5 keV, in agreement with 
that measured above. Due to the hard PSR spectrum, assumed 
to extend up ~ lOOkeV with the same spectral index, the min- 
imum PWN spectrum at higher energies is softer than the total 
emission, but statistically consistent with the spectrum mea- 
sured by IBIS/ISGRI. 

2.5. Radio 
2.5.1. pulsar 

A search for a radio counterpart to the X-ray pulsar was car- 
ried out using the Parkes 64-m radio telescope of the Australia 
Telescope National Facility, CSIRO Astronomy and Space 
Science, in the 10, 20 and 50 cm bands. Observations were 
made with the centre beam of the 20 cm Multibeam receiver 



and analogue filterbank system dManchester et al. 12001b with 
sampling interval of 250 /-is on 2009 October 11 (2 h) and 
2009 November 26 (4 h), and the 10-50 cm receiver and 
Parkes digital filterbank systems (1 h at each of 10 cm and 
50cm), on 2009 October 13. The received centre frequencies 
were 3100, 1374 and 732 MHz at 10, 20 and 50 cm respec- 
tively, and the corresponding bandwidths were 1024, 256 and 
64 MHz. A search of the October data over a period range of 
+60 /is about the nominal pulsar period (see Table [TJ and a 
dispersion measure (DM) range of to 1500 cirT 3 pc yielded 
some possible candidates with DMs in the range 400 - 600 
crrT 3 pc and signal-to-noise (S/N) ratios of < 6.5. A similar 
analysis of the longer 20 cm observation in November gave a 
positive detection with S/N ratio of 8.8 at the predicted pulse 
period and at a DM of 563 ± 4 cm -3 pc, as shown in Figure[3] 
This detection corresponds to one of the candidates from the 
analysis of the October 20 cm observations. 




-0.1 -0.05 0.05 0.1 

Delta P (microsec) 

Fig. 5. — Plot of the 20 cm radio pulse S/N ratio in the pulse period - DM 
plane around the detection values from the 2009 November Parkes observa- 
tions. The greyscale is linear between S/N values of and 10. S/N contours 
at 5, 6, 7 and 8 are overlaid. 



Figure [6] shows the 20 cm (1374 MHz) pulse profile from 
the November observations and the 10 cm (3100 MHz) pulse 
profile obtained by folding and dedispersing the data with the 
nominal pulsar parameters and a DM of 563 cm" 3 pc. The 
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Pulse Phase 

Fig. 6. — Radio pulse profiles at 10 cm (upper) and 20 cm (lower) folded at 
the predicted pulse period and dedispersed with a DM of 563 cirT 3 pc. The 
pulse has been approximately centered in both plots. 



pulsed emission is clearly detected at 10 cm although the S/N 
ratio is low. The profile appears very broad and almost si- 
nusoidal at both frequencies. The profiles have been cali- 
brated in flux density units, giving mean pulsed flux densi- 
ties of approximately 0.25 mJy and 0.11 mJy at 20 cm and 
10 cm, respectively. The upper limit on the pulsed flux den- 
sity at 50 cm is about 0.5 mJy. At face value, these re- 
sults imply a rather flat pulsar spectrum with spectral index 
a ~ -1.0 (throughout this paper, S v oc v"). Radio observa- 
tions with improved sensitivity and with contemporaneous X- 
ray observations will be required to explore the details of the 
radio pulse morphology and the phase relationship with the 
X-ray pulse. Distance estimates ba sed on the observed DM 
are 9.4 kpc dTaylor & Cordes| [T993l) and 11.3 kpc (NE2001, 
iCordes & Lazioll2002l) . Given the large uncertainties inherent 
in DM-based distances, a reasonable estimate for the pulsar 
distance is 10+3 kpc. 

2.5.2. wind nebula 

Young PSRs are often associated with wind nebulae, usu- 
ally detected in the radio and X-ray domains. We searched for 
the radio counterpart of the Chandra SNR in the archival data 
of the Molonglo Galactic Plane Survey at 843 MHz (MGPS-2, 
Murphy et al. 2007), the 2.4 GHz Parkes su rvey of the South- 
ern Galactic Plane dDuncan et al. 1 1 19951) , and the Parkes- 
MIT-NRAO (PMN) survey at 4 85 GH z dGriffith & Wrightl 
119931: ICondon. Griffith & Wrightl 11993b . We found an ex- 
tended MGPS-2 source, MGPS J140045-632542, located at 
R.A.(J2000) = 14 h 00 m 45.17 s , Decl.(J2000) = -6325'42.2", 
3.5" away from the Chandra point-like source. The intrin- 
sic source size of 44.8 x 31.2" (after deconvolution from the 
MOST beam size of ~ 45") is depicted in Figure Q] (left) by 
the dashed ellipse, which nicely matches the X-ray PWN, and 
thus corresponds to its radio counterpart. The 843 MHz flux 
densi ty was measured to be 217.4 ± 9.4 mJy dMurphv et al. I 
12007b . In the 4.85 GHz PMN data, we found an other cata- 
loged source coincid ent with the Chandra SNR. According to 
IWrightet all d 19941) . PMN J1400-6325 lies at R.A.(J2000) 
= 14 h 00 m 45.2 s , Decl.(J2000) = -6325'43". With an angular 
resolution of 5' (FWHM), the source is unresolved, and fea- 
tures a flux density of 113 (131) ± lOmJy, depending on the 



fit pro cedure used (see lGriffith & Wrightlll993UWright et al. I 
119941) . No signi ficant emission is de tected in the Parkes im- 
age at 2.4 GHz dDuncan et al. |[l995l) . and we calculated a 3 
cr upper limit of 0.6 Jy. A power-law fit on these three ra- 
dio measurements gives a spectral index a of -0.33 + 0.05, 
-0.205 + 0.008, and a flux density at 1 GHz of (0.207-0.204) 
± 0.008 Jy, depending on the adopted MGPS-2 source flux. 
Such a index is similar to tho se measured in other radio 
PWNe dGaensler & Slandl2006T) . and we can safely conclude 
these MGPS-2 and PMN sources are the counterparts of the 
X-ray wind nebula discovered with Chandra. 

No radio shell-type structure is detected in the MGPS-2 im- 
age, even though the MOST angular resolution should allow 
us to separate it from the PWN. We have calculated an upper 
limit on the SNR surface brightness by modelling the SNR as 
a uniform annulus with a mean radius of 1.15' and a width 
of 0.15' (see Fig. Q] right) and convolving it with the MOST 
beam size. The upper limit, derived by summing all the pix- 
els within a region encompassing the FWHM of the simulated 
MGPS-2 shell emission, is 2.5 x 10~ 21 W m -2 Hz -1 sr _1 at 
843 MHz, which corresponds to 2.3 x 10~ 21 W irT 2 Hz 1 sr 
1 at 1 GHz (for a = -0.5), as com monly used in studies of 
radio shell-type SNRs dGreenll2009l) . 

3. DISCUSSION 

3.1. General considerations 

After the Crab pulsar, PSR J 1400-6325 is one of the most 
energetic pulsars known in the Galaxy, rivaling the recently 
discov ered PSRs associated with G12.82-0.02 and GO 9+0. 1 
SNRs dGotthelf & Halpernl 120091: ICamilo et al. 1 12009|) . Its 
characteristic age t c of 12.7 kyr likely overestimates the 
true age t, based on comparison with the SNR estimate 
(see below), consistent with the fact that the pulsar has not 
spun down significantly. We then assume in the follow- 
ing that r «c t c {i.e. P ~ Po and £ ~ Eo). The com- 
monly used relations between PSR/PWN X-ray spectra (2- 
1 keV indices and luminosi t ies) and the intr i nsic PSR prop- 
erties dPossenti et al. I l2002t iGotthelfl 12003b iMattana et al~l 
l2008llLiet al. II2008I) . provide consistent estimates with those 
measured here, though with the large uncertainties inher- 
ent of such correlation studies. Amongst these relations, 
those between PSR/PWN X-ray luminosities and £ point 
towards a distance to the source of 6-10 kpc, compati- 
ble with the DM-based estimate of 10 ± 3 kpc. On one 
hand, the high Nh derived from the Chandra spectra sug- 
gests either a large source distance, or local absorption that 
could be undetectable in the existing large-scale HI and 
12 CO Galactic surveys at moderate angular resolutions. Us- 
ing the latest relation between optical extinction and Nh of 
iGiiver & Ozell d2009l) . we find an extinction Ay ~ 9 toward 
IGR J 14003-6326. Further more, according to the Galactic 
interstellar extinction map of [Marshall e t al. ( 2006), such A v 
along the IGR J14003-6326 line of sight corresponds to a 
distance of 5-8 kpc. For such a high extinction, only in the 
case of an over-luminous core-collapse SN event (with My 
< -17, see Richardson 2009), would the stellar explosion of 
G310.6-1.6 have been visible with naked eye (my < 6) for 
distances > 5 kpc. On the other hand, its Galactic latitude b ~ 
- 1 .6 translates into an height from the Galactic Plane of ~ 
280 pc at 10 kpc, which is much larger than the typi cal scale 
heigh t of the Galactic molecular disk (~ 55 pc, e.g. iFerrierd 
2001) within which core-collapse SNe are known to usu- 
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ally explodeQ. We also note that the Galactic longitude of 
G3 10. 6- 1. 6 correspon ds to the Crux-Scutum arm tangent, at 
6-7 kpc (Vallee 2008), where much of the molecular mate- 
rial along this line of sight lies. In order to account for all 
the above-mentioned estimates with their respective uncer- 
tainties, we then adopt a distance of 7 kpc, which implies a 
height similar to that of the Crab nebula, and scale the dis- 
tance in terms of dj = d/7 kpc. 

The very small angular size of G3 10. 6- 1.6 ~ 2.3' = 4.7 dq 
pc already places it amongst the smallest known SNRs, and 
could in turn imply a very young age, similar to the re - 
cently rediscovered Gl. 9+0.3 SNR (iReynolds et al. 1 120081) . 
and the well known Tycho, Kepler and Cas A SNRs. For a 
uniform ambient medium of density no cirT 3 , the SNR has 
swept up only M sw = (47r/3)«o^ 3 ~ 1-3 nodi M Q and thus 
may still be in free expansion phase. In such case, for a 
typical ejected mass of M e j = \0M\$M Q and a SN explo- 
sion kinetic energy of 10 51 £51 ergs, the expansion velocity 
would be ~ 3200(£ , 5i/Mio) 1/2 kms- 1 , implying an age ~ 720 
dj{M\o/Esi) ' 2 yr. Such a young age is also supported by the 
nearly circular shape of the SNR, with the PSR lying very 
near its geometrical center. An offset of < 2.5", as estimated 
on the Chandra image, implies a projected space PSR veloc- 
ity of only < 80 ^tJ 1 kms" 1 , with t = 10 T 3 yr, i.e. sig- 
nifican tly smaller than the average found by IHobbs et al. I 
( 12005b for young PSRs. If we assume G3 10.6- 1.6 has pro - 
duced the same amount of 44 Ti as Cas A dRenaud et al. 1200*61) . 
the IBIS/ISGRI upper limit on the two low-energy 44 Ti lines 
would translate into a lower limit on its age of ~ 280 yr at 
7 kpc. If younger, G3 10. 6- 1.6 would not have originated 
from one of these 44 Ti-rich SN explosions whose remnants 
(except Cas A) have been unfruitfully searched by several 
missions over t he past decades (e.g. iRenaud et al. I l2004t 
iThe et al. 112006b . 

3.2. PWN broadband emission and energetics 

In order to account for the different spectral ind exes mea- 
sured in radio and X-ray domains (see section |2~31 and Table 
13, the PWN broadband synchrotron spectrum should exhibit 
a break at ~ 3 x 10 4 GHz. Such frequency, if originated from 
synchrotron cooling of accelerated electrons inside the PWN, 
would correspond to an uncomfortably high magnetic field 
strength B ~ 330 t" 2/3 /jG. In such case, the PWN should 
have appeared smaller in X-rays than in radio, in contradic- 
tion with what is shown in FigureQ] This break could instead 
reflec t the injected electron spectrum (see e.g. ICamilo et al. I 
120061) or the nature of th e magnetic field and of the re sult- 
ing synchrotron emi ssion dFleishman & Bietenholzl2007l) . As 
presented in section [2~4l the break measured at ~ 6 keV could 
be that expected from synchrotron cooling and, thus, provides 
a PWN age t ~ 880 B~^ 2 yr, with B = 10B 10J uG. VHE obser- 
vations can help further constrain t he magnetic fi e ld stre ngth 
inside the PWN. As reported by IChaves et all ( 12008b . no 
point-like source toward IGR J 14003 -6326 is detected above 
5 cr in the latest map of the H.E.S.S. Galactic Plane Survey, in 
15 hours of effective time. This tra nslates into a H.E.S .S. up- 
per limit of 4 % of the Crab nebula dKhelifi et al. 12008b . i.e. 3 
x 10~ 12 ergs crrT 2 s _1 in the 1-10 TeV energy range. In a 
simple one-zone leptonic scenario, the VHE emission comes 

8 We note that G310.6-1 .6 could po tentially originate from a runaway OB 
/ Wolf-Rayet star (see e.g. IDra v et al. 2005). This scenario would explain 
the large height from the Galactic Plane for a distance of 10 kpc. 



from inverse-Compton scattering of VHE electrons on ambi- 
ent photons (CMB and infrared and optical Galactic fields, 
from dust and stars, respectively). From the ratio between the 
X-ray flux and the H.E.S.S. upper limit, we find that the mag- 
netic field must be larger than 6 /jG, and then, the age must be 
smaller than ~ 1900 yr. The corresponding maximal energy 
content in particles E p max amounts to ~ 2.2 x 10 48 d 2 , ergs. 

PWN energetics can also provide useful insight on the 
source age, as young PWNe usually feature a ratio between 
the internal energy E mt and the PSR rotational spin-down en- 
ergy E sc i = fir < 1 (~ 0.44 for r <K the initial PSR spin-down 
time To, and depends w eakly on the supernova density pro- 
file, see IChevaheril2"00l equation 30). On one hand, E n max 
thus provides a conservative limit E s ,; mav ~ 5 x 10 48 c/ 2 ergs, 

which implies a maximal age of ~ 3000 yr for E = 5.1 x 10 37 
ergs s _1 . On the other hand, a lower limit on E mt can be ob- 
tained from equipartition arguments dGovoni & Ferettill2004t 
IChevalierJl2005T) . We estimate that £i n t,min ~ 4 x 10 46 d' 7/7 
ergs, i.e. «: £t3 = 1.6 x 10 48 ergs, and a magnetic field at 
equipartition B eq ~ 81 af" 4 ^ 14 /iG, from the radio properties 
of the PWN. In case IGR J 14003-6326 is at the equipartition 
between particles and magnetic field, its age would then be 
<§: 1 kyr. However, H.E.S.S. has recently discovered VHE y- 
ray emis sion toward two other you ng PWNe, G2 1.5 -0.9 and 
Kes 75 (iDiannati-Atai et aT7ll2008h . These observations im- 
ply magnetic fields of 10-15 fiG, well below the eq uipartition 
value s, and therefore support earlier suggestions ( IChevalierl 
12001 that some young PWNe are particle dominated. If 
the same applies to IGR J 14003-6326, E[ nt at the equipar- 
tition does not reflect the whole PWN energy content, and 
£int(B) °c B~ 3/2 in the particle-dominated regime. In this 
regime, following IChevaUed (12005b . the E m (d,B)/ET ~ 0.44 
relation, valid when t <k to as assumed here, would be con- 
sistent with the existence of a break measured at ~ 6 keV (t 
~ 880 B~q 2 yr) for a distance of 7-8 kpc. 

3.3. The evolutionary stage of G310.6—1.6 

As discussed previously, several arguments point towards a 
very young composite SNR. The ratioj between the PWN and 
SNR radii is R « 1/4, and can be used to probe the parame- 
ters governing th e early evolution of composite SNRs, namely 
E51, M e j and «o (Ivan der Swaluw et al. Il200lt iBlondin et al. I 
l200l . given the £ measured with RXTE. First, «o and £51 
are constrained by the faint X-ray level, and the 1 GHz up- 
per limit, of the SNR synchrotron surface brightness. Ac- 
cording to lBerezhko & Voikl (120041) . such low values can only 
be explained by a very tenuous medium (no < 0.01 crrT 3 
together with a sub-energetic explosion. Furthermore, usin g 
the so-called radius method of lvan der Swaluw & Wul J200l . 
the rotational energy £«/ relates to the PWN to SNR radius 
ratio R as: E sd ~ 2 x 10 48 £51 (R/0.25) 3 (t] 3 /2)- 3 erg, or 
t ~ 1200 £51 (R/0.25) 3 (773 /2) -3 yr, with 773 a dimensionless 
parameter ranging between 1 and 3 during the early evolution 
phase (see their Fig. 1). With the upper limit on E st / estimated 
from the H.E.S.S. non detection, and the lower limit £i n t, m in 
from equipartition arguments, we infer 0.02 dV^ (773 /2) 3 < 

9 The radio PWN found in the MGPS-2 data is elliptical, so that we sup- 
pose its effective size to be V44.8 X 31.2". This yields a PWN to SNR size 
ratio of ~ 37"/2.3' ~ 1/4. 

10 This argues against a local absorption at the origin of the large Nh, and 
distances much less than the Crux-Scutum spiral arm tangent. 
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£51 < 1 (773 /2) 3 , in agreement with the weak radio and X- 
ray surface brightnesses of the SNR. We estimate the initial 
PSR period Po > 29.5 ms (for a braking ind ex n - 3), for 
a curr ent p eriod of 3 1 18 m s. Based on the iBlondin et al. I 
(120011) and iChevalierl d2005l) calculations, we then explored 
a wide range of £51, no and M e j (the power-law index of the 
ejecta distribution was set to 12). For a distance of 7 kpc, if 
no is set to 0.01 crrT 3 , we find that only a sub-energetic (£51 
= 0.05), low ejecta mass (M e j = 3M ) SN explosion could 
explain the measured R. The resulting SNR age and PWN 
magnetic field are ~ 920 yr and ~ 10 //G, respectively. Note 
that the constraint on «o is only qualitative, as a quantitative 
estimate would require a detailed modelling of the SNR syn - 
chrotron emission, as performed bv lBerezhko & Volkl ( 120041) . 
If we relax this constraint and set no to the Galactic value of 
1 cm -3 , this would result in a larger energy of the explosion 
(£51 — 0.4), and a smaller SNR age (t ~ 550 yr), for the same 
mass of the ejecta. In this regard, knowing the nature of the 
SNR X-ray emission is crucial to assess the thermal flux level, 
and then, the surrounding density. 

4. CONCLUSION 

Even though some assumptions have been made in the 
previous estimates (e.g. constant PWN magnetic field and 
PSR spin-down power), they give valuable insight on the na- 
ture of this new Galactic composite SNR. G3 10. 6- 1.6 har- 
bors a highly energetic 31.18 ms pulsar, PSR J 1400-6325, 
which powers a wind nebula, both lying at the center of 
the host shell. All of the existing multi-wavelength obser- 
vations suggest it is a young SNR (< 10 3 yr), and most 
likely distant (> 5 kpc). However, many questions still re- 
main to be answered. First, the distance is not very well 
constrained, though a large distance is favored by the dis- 
persion measure of the radio pulse emission and by the 
large Nh measured with Chandra. High-resolution observa- 
tions of the ISM tracers such as HI and 12 CO are then war- 
ranted to assess the surrounding medium properties. More- 
over, the break measured at 6 keV and the spectral soft- 
ening at increasing distances in the PWN need to be con- 
firmed, and the SNR X-ray spectrum needs to be investi- 
gated with more X-ray data. Nevertheless, G310.6-1.6 falls 
into the emerging class of "multi-wavelength", young (t < 
a few 10 3 yr) and composite SNRs, harboring very en- 



ergetic PSRs and wind nebulae shining in radio, X- 
rays and potentially in HE/VHE y-rays. The list includes 
Kes 75 and G2 1.5-0.9 (IGotthelf et al. I HOQOt ICamilo et al. I 
l2"006t iBietenholz & Bartell 120081: iDiannati-Atai et al. 112008b . 
and more recently G0.9+0~l (lAharonian et al. I I2005al) and 
HESS J1813-178 (lAharonian et al. I I2005b| ) ~whose long- 
expec t ed PSRs have recen t ly bee n discovered (ICamilo et al. I 
120091; IGotthelf & Halpernl 120091) . It is of inte rest to note 
that Fe rmi/LA^T has n ot detected a bright source (lAbdo et al. I 
l2009ah . nor pulsar ( lAbdo et al. I l2009bllcT) coincident with 
G310.6-1.6, although, at first glance, energetic PSRs should 
be the most easily detectable sources. However, not all of 
the above-mentioned young and energetic PSRs have been 
detected by Fermi. From the first Fermi catalog of gamma- 
ray PSRs, the sensitivity for a blind search of pulsed emission 
in the Galactic Plane is co nservatively taken t o be 2 x 10~ 7 
cirT 2 s _1 above 100 MeV dAbdo et al H2009d) . Assuming a 
spectrum similar to that of PSR J1833— 1034 associated with 
G21. 5-0.9, PSR J1400-6325 features a maximal efficiency 



£ r /E of ~ 1 di. 



This is close to what is measured 



from other young PSRs (see Fig. 6 of lAbdo et al~l d2009d) ). 
and argues in favor of a fairly large distance to the source, 
as outlined in section [3] Further HE/VHE observations of 
G3 10. 6- 1.6 will certainly provide important constraints both 
on the PSR gamma-ray spectrum and on the PWN magnetic 
field strength. 
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